NUWC-NPT  Technical  Report  10,632 
27  January  1997 


A  Comparison  of  Vector-Sensing  and 
Scalar-Sensing  Linear  Arrays 

Benjamin  A.  Cray 
Submarine  Sonar  Department 

Albert  H.  Nuttali 

Surfoce  Undersea  Warfare  Directorate 


Naval  Undersea  Warfare  Center  Division 
Newport,  Rhode  Island 

Approved  for  public  release;  distribution  is  unlimited. 

199/0328  056 


PREFACE 


This  report  was  prepared  xmder  the  “Submarine  Large  Passive 
Flank  Array  Project”  sponsored  by  the  Office  of  Naval  Research  (ONR); 
Technology  Project  RJ14A15.  The  NUWC  Division  Newport  Project  No. 
is  H60030;  NUWC  project  manager,  Howard  H.  Schloemer  (Code  2002); 
principal  investigator,  Benjamin  A.  Cray  (Code  2133).  The  ONR  Project 
Manager  is  Robert  Varley,  Tactical  Systems  Science  and  Technology 
Program  (Code  321SS). 

Some  of  the  work  described  here  was  co-sponsored  by  the 
Independent  Research  (IR)  Program  of  NUWC,  Project  B 10007,  entitled 
“Near-Optimum  Detection  of  Random  Signals  with  Unknown  Locations, 
Structure,  Extent,  and  Strengths.”  The  IR  program  is  funded  by  the  Office 
of  Naval  Research;  the  NUWC  Division  Newport  program  manager  is 
Stuart  C.  Dickinson  (Code  102).  The  principal  investigator  for  this  IR 
project  is  Albert  H.  Nuttall  (Code  311). 

The  technical  reviewer  for  this  report  was  Ronald  R.  Kneipfer 
(Code  214). 


Reviewed  and  Approved:  27  January  1997 


R.  J.  Martin 

Head,  Submarine  Sonar  Department 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  respi 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  infom 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headqua 
Davis  Highway.  Suite  1204.  /Vlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budc 

□nse,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
nation.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
rters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
et.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  Blank) 

2.  REPORT  DATE 

27  January  1997 

3.  REPORT  TYPE  AND  DATES  COVERED 

Final 

4.  TITLE  AND  SUBTITLE 

A  Comparison  of  Vector-Sensing  and  Scalar-Sensing  Linear  Arrays 

5.  FUNDING  NUMBERS 

PR  H60030 

PR  B10007 

6.  AUTHOR(3) 

Benjamin  A.  Cray  and  Albert  H.  Nuttall 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Undersea  Warfare  Center  Division 

1176  Howell  Street 

Newport,  Rhode  Island  02841 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TR  10,632 

9.  SPONSORING/MONITORING  AGENCY  NAME{S)  AND  ADDRESS(ES) 

Office  of  Naval  Research  (Code  321SS) 

800  North  Quincy  Street 

Arlington,  Virginia  22217 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

Some  of  the  features  which  are  unique  to  beamforming  a  linear  array  of  vector  sensors,  not  present  with  scalar-sensing 
elements  (such  as  conventional  pressure  sensors),  are  described  in  this  report.  Vector  sensors  measure  acoustic  particle 
motion,  for  example  acoustic  particle  velocity,  in  a  given  direction  and  these  vectors  have  complex-valued  components  which 
are  proportional  to  the  amplitude  and  phase  of  a  propagating  acoustic  field.  Three  types  of  vector  sensors  are  considered;  a 
uniaxial  vector  sensor  which  measures  acoustic  particle  velocity  along  a  single  axis,  a  biaxial  vector  sensor  measuring  velocity 
in  two  orthogonal  directions,  and  a  triaxial  sensor  which  measures  all  three  orthogonal  components  of  the  velocity  vector. 

Comparisons  are  made  between  the  performance  (in  terms  of  received  element  and  array  sensitivity,  array  directivity, 
and  signal  gain  with  array  curvature)  of  arrays  comprised  of  each  type  of  vector  sensor,  to  arrays  comprised  of  pressure 
sensors  only.  The  comparisons  have  also  evaluated  two  boundary  conditions,  perfectly  baffled  sensors  and  sensors  in  the  free 
field.  It  is  shown  that  a  single  triaxial  vector  sensor  may  be  steered,  in  both  azimuth  and  elevation,  to  any  planewave  arrival 
direction.  A  single  pressure  sensor  provides  no  directivity  in  the  free  field.  A  comparison  of  array  directivity  indicates  (for 
optimal  real-valued  amplitude  shading  coefficients)  that  triaxial  and  biaxial  vector  sensors  can  provide  5  dB  of  directivity  gain 
over  that  obtained  from  an  identical  array  of  pressure  sensors.  The  gain  varies  with  the  direction  the  array  is  steered  to. 
Similarly,  multiaxis  vector  sensors  improve  a  conformal  array's  signal  gain  for  all  array  steerings.  For  a  cyiindrically-curved 
vertical  array,  vector  sensors  improved  signal  gain  by  a  minimum  of  approximately  3  dB  at  45  degree  steerings,  to  a  maximum 
gain  of  4  dB  at  broadside. 

It  has  also  been  shown  that,  unlike  pressure-sensors,  there  exist  no  interelement  array  spacing  at  which  vector  components 
are  all  uncorrelated  with  themselves  for  a  spherically  isotropic  field.  Thus,  even  for  half-wavelength  spacing,  element  noise 
correlations  must  be  accounted  for  in  any  signal  detection  procedure. 


14.  SUBJECT  TERMS 

Acoustic  Vector  Sensors 

Submarine  Sonar  Systems 

15.  NUMBER  OF  PAGES 

40 

Scalar  Pressure  Sensors 

Velocity  Sensors 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  OF  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  ABSTRACT 

Unclassified  Unclassified  Unclassified  SAR 


TABLE  OF  CONTENTS 


Page 


LIST  OF  ILLUSTRATIONS .  ii 

INTRODUCTION .  1 

REVIEW  OF  PREVIOUS  RESEARCH .  1 

MODEL  OF  ARRAY  RESPONSE . 2 

Amplitude  Response  of  Vector  Sensors  .  5 

Beamforming  a  Single  Vector  Sensor .  10 

Beamforming  an  Array  of  Vector  Sensors . 12 

NUMERICAL  RESULTS  .  16 

Comparison  of  Beam  Patterns  .  16 

Free-Field  Directivity  Calculations .  18 

The  Effect  of  Array  Curvature .  22 

CONCLUSIONS .  24 

APPENDIX  A  —  ANALYTICAL  EVALUATION  OF  DIRECTIVITY  FOR 

HALF- WAVELENGTH  ARRAY  SPACING .  A-1 

APPENDIX  B  —  DERIVATION  OF  MAXIMUM  ARRAY  GAIN .  B-1 

APPENDIX  C  —  SENSOR  CORRELATIONS  IN  AN  ISOTROPIC  NOISE  FIELD C-1 

REFERENCES  .  R-1 


1 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Array  Coordinates  and  Geometry  for  an  Equispaced  Line  Array  of 

Multiaxis  Vector  Sensors .  3 

2  Angular  Coordinates  Used  to  Define  Incident  Direction  of  Acoustic 

Planewave  at  Wavevector  k  .  5 

3  Triaxial  Vector  Sensor  Element  Response  on  an  Ideal  Boundary 

Showing  Conical  Angles  0^,  ©y,  0^ .  7 

4  Amplitude  Response  of  a  Half-Omnidirectional  Pressure  Sensor 

on  a  Rigid  Boundary .  7 

5  Amplitude  Response  of  the  v-^-Component  of  a  Uniaxial  Vector 

Sensor  on  an  Ideal  Pressure  Release  Boundary  .  7 

6  Cross  Section  of  the  Amplitude  Response  of  the  Vy-Component  of  a 

Uniaxial  Vector  Sensor  versus  Azimuth  Angle  (((>  =  0°) .  8 

7  Amplitude  Response  of  the  V;c-Component  of  a  Uniaxial  Vector 

Sensor  on  an  Ideal  Pressure  Release  Boundary . . .  9 

8  Amplitude  Response  of  the  v^-Component  of  a  Uniaxial  Vector 

Sensor  on  an  Ideal  Pressure  Release  Boxmdary .  9 

9  Magnitude  of  the  Vjc-  and  v-y-Components  of  a  Biaxial  Vector 

Sensor  on  an  Ideal  Pressure  Release  Boundary .  9 

10  Magnitude  of  the  y^;-  and  y^-Components  of  a  Biaxial  Vector 

Sensor  on  an  Ideal  Pressure  Release  Boundary  .  9 

1 1  Magnitude  of  the  v^-  and  y^-Components  of  a  Biaxial  Vector  Sensor 

on  an  Ideal  Pressure  Release  Boundary .  10 

12  Magnitude  of  the  Vx-,  vy,  and  y^-Components  of  a  Triaxial  Vector  Sensor 

on  an  Ideal  Pressure  Release  Boimdary .  10 

13  Amplitude  Pattern  of  a  Single  Triaxial  Vector  Sensor  Steered  to 

45°  Azimuth  and  30°  Elevation  .  1 1 

14  Beam  Pattern  of  a  Single  Triaxial  Vector  Sensor  Steered  to 

45°  Azimuth  and  30°  Elevation  .  11 

15  Beam  Pattern  in  Elevation  of  a  Single  Triaxial  Vector  Sensor  Steered 

to  Elevation  Angles  of  0°,  45°,  and  90° .  12 


11 


LIST  OF  ILLUSTRATIONS  (Cont’d) 

Figure  Page 

16  Beam  Pattern  at  Broadside  of  a  Linear  10-Element  Triaxial 

Vector-Sensing  Array  Having  Uniform  Shading  Weights  .  16 

17  Beam  Pattern  at  Broadside  of  a  Linear  10-Element  Pressure- 

Sensing  Array  Having  Uniform  Shading  Weights .  16 

18  Comparison  of  the  Azimuthal  Beam  Pattern  of  Vector-  and 

Pressure-Sensing  Arrays  at  Broadside .  17 

19  Beam  Pattern  of  a  Linear  10-Element  Triaxial  Vector-Sensing  Array  Having 

Uniform  Shading  Weights  Steered  to  45°  Azimuth  and  30°  Elevation .  17 

20  Beam  Pattern  of  a  Linear  10-Element  Pressure-Sensing  Array  Having 

Uniform  Shading  Weights  Steered  to  45°  Azimuth  and  30°  Elevation .  17 

21  Maximum  Directivity  for  an  N-Element  Equispaced  Line  Array  at 

Broadside  (Oj  =  0°)  Utilizing  Optimal  Real  Weights .  19 

22  Maximum  Directivity  for  an  N-Element  Equispaced  Line  Array  at 

Endfire  (Gj  =  90°)  Assuming  Optimal  Real  Weights .  19 

23  Directivity  Versus  Azimuthal  Angle  for  a  10-Element  Equispaced 

Line  Array  at  Elevation  Angle  ((t>s  =  0°)  .  20 

24  Directivity  Versus  Azimuthal  Angle  (Gj)  for  a  10-Element  Equispaced 

Line  Array  at  Elevation  Angle  ((j)3  =  15°)  .  21 

25  Directivity  Versus  Azimuthal  Angle  (G^)  for  a  10-Element  Equispaced 

Line  Array  at  Elevation  Angle  (<|)5  =  30°)  .  21 

26  Section  of  a  Curved  Vertical  Line  Array  Comprised  of  Biaxial  Vector 

Sensors  (v^,  v^)  and  Component  Normals  («j,',«z') .  22 

27  Signal  Gain  of  a  Fully-Populated  Cylindrically  Curved  Vertical  Line  Array  of 

Biaxial  Vector  Sensors  with  Uniform  Shading  Weights  .  24 


iii/iv 

Reverse  Blank 


A  COMPARISON  OF  VECTOR-SENSING  AND 
SCALAR-SENSING  LINEAR  ARRAYS 


INTRODUCTION 

Unlike  scalar  pressure  sensors,  acoustic  vector  sensors  measure  the  amplitude  and  phase 
of  acoustic  particle  motion  in  a  given  direction.  These  sensors,  generally  referred  to  as  velocity 
sensors,  measure  a  vector  quantity,  i.e.,  acoustic  particle  acceleration,  velocity,  or  displacement. 
A  single-axis,  or  uniaxis,  acoustic  vector  sensor  will  be  defined  here  as  a  sensor  that  measures 
one  (cartesian)  component  of  the  acoustic  field  vector.  Biaxial  vector  sensors,  such  as  biaxial 
accelerometers,  denote  sensors  that  measure  motion  in  two  orthogonal  directions,  while  triaxial 
vector  sensors  measure  all  three  orthogonal  components.  These  uniaxis  and  multiaxis  vector 
sensors  can  then  be  configured  into  an  array  of  elements,  the  elements  of  which  may  be  delayed, 
weighted,  and  summed,  as  commonly  done  with  standard  pressure  hydrophones.  However,  Aere 
are  some  unique  features  present  in  Ae  beamforming  of  vector  sensors  that  are  not  present  with 
pressure  sensors.  This  report  will  describe  some  of  these  unique  features,  i.e.,  the  differences 
between  arrays  of  scalar  pressure  sensors  and  arrays  composed  of  acoustic  vector  sensors. 


REVIEW  OF  PREVIOUS  RESEARCH 

A  brief  literary  review  indicates  that  a  group  of  researchers  have  recently  documented 
some  aspects  of  beamforming  arrays  of  multiple-axis  vector  sensors.  D'Spain^  constructed  and 
deployed  a  16-element  vertical  line  array  that  measured  the  three  components  of  fi:ee-field 
acoustic  particle  velocity  as  well  as  acoustic  pressure.  The  low  firequency  Directional  low 
Frequency  Analysis  and  Recording  (DIFAR)  array  consisted  of  collocated  pressure  hydrophones 
with  three  orthogonally  mounted  geophones.  The  output  firom  a  geophone  sensor  is  directly 
proportional  to  acoustic  particle  velocity.  At-sea  data  collected  fi'om  the  array  were  beamformed 
conventionally  (delayed  and  summed  with  uniform  amplitude  weighting)  and  adaptively,  using  a 
minimum  variance  Capon  approach  on  all  four  quantities. 

Nehorai^  developed  an  analytical  model  to  compare  the  direction-of-arrival  (DO A) 
estimation  performance  of  an  array  of  sensors,  which  measures  acoustic  pressure  and  all  three 
components  of  acoustic  particle  velocity,  to  that  of  an  array  which  measures  acoustic  pressure 
only.  The  paper  treated  both  single  and  multiple  acoustic  sources  in  an  ideal  fi-ee-field 
environment.  Cramer-Rao  lower  bounds  (CRLB)  on  DOA  estimation  for  the  single-source 
single- vector  sensor  as  well  as  multiple-source  multivector  sensors  were  given.  The  CRLB 
provides  a  universal  measure  of  achievable  estimation  accuracy  and  is  commonly  used  by  signal 
processors  to  compare  the  merits  of  different  high-resolution  DOA  methods.  Hawkes^  extended 
the  work  of  Nehorai  by  introducing  a  plane  boundary  at  or  near  the  vector-sensing  array.  As 
done  previously  in  the  paper  by  Nehorai,  expressions  for  the  CRLB  were  derived  and  analyzed 
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for  various  stand-off  distances  from  ideal  pressure-release  and  rigid  boundaries.  In  reference  4, 
Kneipfer,  using  notation  and  formulations  common  to  the  sonar  community,  rederived  some  of 
the  analytical  model  developed  by  Nehorai.  The  work  reported  here  is  a  continuation  of 
Kneipfer's  preliminary  efforts,  examining  in  additional  detail  the  merits  of  vector-sensing  arrays. 


MODEL  OF  ARRAY  RESPONSE 

This  section  will  develop  a  model  for  the  beamformed  response  of  uniaxial,  biaxial, 
and  triaxial  arrays  of  vector  sensors  which  measure  acoustic  particle  velocity  in  the  free  field  and 
on  an  idealized  boundary.  The  boundary  is  assumed  to  provide  perfect  baffling  of  half  the  free- 
field  acoustic  field  and  ideal  pressure-release  behavior  for  each  orthogonal  acoustic  velocity 
component.  The  baffled  and  free-field  vector  array  responses  are  compared  to  that  of  a  similar 
array  of  scalar  pressure  sensors.  The  effects  of  a  nonideal  boundary,  such  as  a  submarine  hull 
coating,  may  be  incorporated  in  future  investigations.  It  is  important  to  stress  that  these 
comparisons  have  been  made  in  order  to  assess  the  merits  of  an  array  of  vector  sensors  to  an 
array  of  scalar  sensors  based  on  ideal  conditions.  If  the  benefits  of  vector-sensing  are  significant, 
then  subsequent  investigations  can  incorporate  realistic  boundary  behavior.  Keltie^  has 
investigated  a  realistic  solid-layer  model  for  typical  submarine  hull  coatings,  which  models  both 
the  coating's  tangential  and  normal  signal-induced  velocities.  If  warranted.  Ids  results  may  be 
incorporated  into  the  model  developed  here. 

Each  multiaxis-vector  sensor  will  be  assumed  to  be  collocated,  measuring  three 
orthogonal  components  of  particle  velocity  at  the  same  spatial  location.  Such  a  restriction  may 
not  be  optimal  for  detection  but  is  used  here  to  simplify  Ae  formulation.  Figure  1  illustrates  the 
array  geometry  and  coordinate  system. 

Acoustic  particle  velocities  are  defined  as  small  amplitude  motions  of  constant  volume 
fluid  particles  (each  particle  having  unvarying  fluid  properties)  about  some  equilibrium  position. 
Acoustic  pressure  is  defined  as  the  instantaneous  pressure  at  a  point,  minus  the  hydrostatic 
pressure  at  that  point.  In  cartesian  coordinates,  time-harmonic  acoustic  planewaves,  which  have 
wavefronts  propagating  in  the  direction  of  acoustic  wavevector  -  k  (by  convention  the  acoustic 
wavenumber  has  been  defined  to  point  from  the  origin  towards  the  source),  can  be  characterized 
by  acoustic  pressure  or  by  acoustic  particle  velocity.  Thus, 


and 


(1) 

(2) 

(3) 

(4) 

where  ||A:||(=  co  /  c)  is  the  acoustic  wavenumber,  with  co  the  circular  frequency,  c  is  the  sound 
speed,  f  =  [x,y,z]  is  the  position  vector,  P  is  the  amplitude  of  the  sound  wave  pressure,  and 
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^5  ^7  are  the  amplitudes  of  the  component  velocities  in  the  x,  y,  and  z  directions, 
respectively.  The  complete  vector  field  may  be  written  as 

v^^^(f,co)  =  v^(r,(o)fi^  +  Vy{r,&)ny  +  (5) 

and  it  is  noted  that  the  velocity  vector  comprises  all  three  vector  components  and  is 
complex;  i.e.,  it  has  complex-valued  components  [Vx,^y,v^-  The  vectors  {n^,ny,ri^}  are 
mutually  perpendicular  unit  vectors  aligned  as  shown  in  figure  1.  The  elements  of  an  uniaxis 
vector-sensing  array  measure  only  one  of  the  velocity  vector  components,  such  as  ,'Vy,  or  v^, 
while  the  multiaxis  arrays  measure  two  or  more  of  these  orthogonal  vectors.  For  example, 

v^^^(f,©)  =  Vj,(f,co)n3„  (6) 

and 

v^^^(f,G))  =  v^(f,co)n^  +  v^{r,(Si)n^.  (7) 

The  superscript  denotes  the  dimension  of  the  velocity  vector,  biaxial  vector  sensors  having  a 
dimension  of  2,  triaxial  a  dimension  of  3.  Other  combinations  are,  of  course,  possible  and  will 
be  discussed  shortly. 


IDEAL 

BOUNDARY 


Figure  1.  Array  Coordinates  and  Geometry  for  an  Equispaced 
Line  Array  of  Multiaxis  Vector  Sensors 
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It  should  be  noted  that,  with  collocated  multiaxis  vector  sensors,  the  component  phases, 
equations  (2),  (3),  and  (4),  may  be  referenced  or  calibrated  in  any  manner  relative 
to  each  other.  That  is,  a  single  multiaxis- vector  sensor  could  be  used  to  determine  whether 
propagating  waves  had  clockwise  or  coimterclockwise  particle  rotations  by  examining  the 
differences  in  the  component  phases.  A  90°-phase  shift  of  a  tangential  component  relative  to  the 
normal  vector  component  would  be  indicative  of  a  surface  wave  having  clockwise  particle 
rotations.  This  is  an  aside,  but  the  relationship  of  the  component  phases  may  prove  to  be  a  useful 
property  in  subsequent  investigations,  such  as  in  the  use  of  vector  sensors  to  cancel  noise  rather 
than  detect  signals.  However,  here  the  components  of  the  vector  sensors  will  be  assumed  to  be 
in-phase,  with  the  component  phases  set  to  zero.  The  harmonic  time-dependence  exp(zcor)  of  the 
acoustic  pressure  and  particle  velocities  will  be  suppressed  henceforth. 

The  spatial  gradient  of  acoustic  pressure  and  the  particle  velocities  are  related  via  Euler's 
linearized  momentum  equation.  Namely,  for  harmonic  waves. 


yp(f)  =  -p^-^  =  z(DPo(v^(r)n^+v^(r)«^+v^(fK),  (g 

where  the  density  of  the  fluid  medium  is  denoted  as  . 

Now,  if  one  considers  a  first-order  Taylor  series  expansion  of  acoustic  pressure  about 
some  point  defined  by  the  position  vector  ,  then. 


p{r)  =  p{r^  )  +  )  •  (r  -  rj 


(9) 


where  equation  (8)  has  been  substituted  for  the  gradient  in  equation  (9)  and  f-r^  =  {x-x^)n^ 

+  (.V  -  Jo )fiy  +{z-z^)n^.  D'Spain  remarked  on  this  expansion,  noting  that  the  simultaneous 
measurement  of  acoustic  pressure  and  particle  velocity  at  a  single  point  is  equivalent  to 
measuring  the  acoustic  pressure  in  a  volume  of  space  about  the  measurement  point.  Therefore,  a 
single  multiaxis-vector  sensor  may  be  considered  similar  to  a  small  (in  terms  of  an  acoustic 
wavelength)  volumetric  array  of  scalar  pressure  sensors,  and  a  line  array  of  multiaxis- vector 
sensors  is  equivalent  to  a  line  array  of  discrete  volumetric  pressure-sensing  elements.  However, 
the  accuracy  of  the  above  expansion,  since  it  is  truncated  to  first  order,  decreases  with  increasing 
volumetric  radius,  R  =  ||f-foll .  This  implies  that  the  ability  of  equation  (9)  to  correctly  measure 
the  acoustic  pressure  in  a  given  volume  of  space  is  dependent  on  the  acoustic  wavelength. 
Comparing  the  ratio  of  the  second-order  Taylor  series  term  (not  shown)  to  the  first-order  term 
given  in  equation  (9)  indicates  that  the  truncated  expansion  is  valid  for  a  volumetric  radius  less 
than  approximately  one-third  of  an  acoustic  wavelength. 

Pressure-gradient  hydrophones  can  measure  the  quantity  Vp  as  well  by  measuring  the 
phase  difference  of  a  planewave  between  two  spatially  separated  hydrophones.  The  separation 
distance  must  increase,  in  order  to  obtain  an  accurate  estimate  of  the  gradient,  with  increasing 
wavelength.  At  low  frequencies,  the  separation  distances  required  for  pressure-gradient 
hydrophones  are  often  too  large  to  be  feasible  for  array  applications. 
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AMPLITUDE  RESPONSE  OF  VECTOR  SENSORS 


The  assumed  geometry  of  an  equispaced  linear  array  of  vector  sensors  is  shown  in  figure 
1.  For  the  baffle  configuration,  the  sensors  lie  in  the  xz-plane  and  are  bounded  by  a  semi-infinite 
acoustic  medium  for  y  >  0.  The  baffle  here  is  assumed  ideal;  i.e.,  perfectly  rigid  in  the  case  of  an 
array  of  pressure  sensors  and  entirely  pressure  release,  for  each  component,  in  the  case  of  an 
array  of  vector  sensors.  In  equation  (1),  P  is  defined  as  the  amplitude  or  angular  sensitivity  of 
the  pressure  sensor  for  angular  coordinates,  (0,<t»),  where  0  is  defined  as  the  azimuthal  angle  and 
(j)  the  elevation  angle.  For  half-omnidirectional  element  sensitivity  on  a  rigid  baffle 


P(0,(t.)  = 


/or-9O°<0,(t)<9O° 

otherwise 


(10) 


It  is  assumed  that  each  orthogonal  component  of  a  vector  sensor  has  an  amplitude 
response  or  element  sensitivity  which  is  proportional  to  the  cosine  of  the  angle  between  the 
planewave  incidence  arrival  direction  and  each  component  normal  (the  component's  mavimnm 
response  axis).  This  assumption  is  based  upon  elementary  solutions  to  the  wave  equation  at  the 
boundary  of  an  ideal  pressure  release  surface.  The  component  normals  are,  for  the  simple  line 
array  configuration  under  investigation,  equivalent  to  the  unit  normals,  {h^ ,  } ,  as  defined 

in  figure  1. 

The  wavevector_describing  the  orientation  of  the  arrival  of  the  incident  planewave,  shown 
in  figure  2,  is  denoted  k  with  components  [kx,  ky,  k^].  Figure  2  also  depicts  the  conical  angle  0y, 
as  in  equation  (16),  left-hand  side,  formed  between  y-unit  normal,  Ry ,  and  the  arrival  direction. 


Figure  2.  Angular  Coordinates  Used  to  Define  Incident  Direction  of 
Acoustic  Planewave  at  Wavevector  k 
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From  figure  1  then, 


=  Wsin(0)cos((|)),  (11) 

=  Wcos(0)cos((t)),  (12) 

and 

=  llBsin((j)).  (13) 

The  angle  formed  between  the  vectors  k  and  {fix  ,ny,nz}  can  be  determined  via  their  dot 
product.  Hence,  cosine  angular  dependence  yields 

T,  .  p  cos(©  • )  for  - 90° < 0,  (})< 90°  and  j =x,y,z 

ji  otherwise 

where 

cos(©^  )  =  cos(<|))sin(0),  (15) 

cos(©  ^ )  =  cos(<t))  cos(0),  (16) 

and 

cos(©2  )  =  sin((t>).  (17) 

The  factor  of  2  in  equation  (14)  is  due  to  the  assumption  of  an  ideal  surface,  thereby 
doubling  the  signal  amplitude  for  all  angles  of  incidence.  Thus,  a  baffled  triaxial  vector  sensor 
can  be  expressed  as 

v^^^(0,<j>,f^)  =  2{cos(©jj)e'**^'’n;^  +cos(©^)e**^^°«^  +cos(©2)e'^^^"«2}’  (18) 

with  the  velocity  vector  dependence  on  azimuthal  and  elevation  angles  now  explicitly  indicated. 

Another  illustration  of  the  conical  angles  associated  with  each  vector  component  is  given 
in  figure  3.  Here  the  amplitude  response  of  the  vector  components,  Vx,  Vy,  and  v^,  on  an  ideal 
boundary,  are  characterized.  In  the  free  field,  of  course,  the  vector  sensor  responses  are  simply 
reduced  by  one-half  and  extended  over  all  acoustic  space  (-180°  <  0  <  180°,  —90°  <^<  90°). 

The  amplitude  response  of  a  half-omnidirectional  pressure  sensor,  given  by  equation  (10), 
is  shown  as  a  gray  scale  image  in  figure  4.  In  all  of  these  images,  the  darker  regions  correspond 
to  minima  and  bright  to  maxima,  and  for  figures  4  through  12,  the  gray  scale  varies  from  zero  to 
a  maximum  value  of  2,  shown  as  white.  The  abscissa  is  azimuthal  angle  0  from  -90°  to  -1-90°;  the 
elevation  angle  is  shown  on  the  ordinate  over  the  same  angular  region.  For  illustration,  in  figure 
4  only,  the  half-omnidirectional  critical  angle  was  set  to  +89°.  As  indicated,  the  amplitude 
response  of  a  half-omnidirectional  pressure  sensor  is  essentially  uniform  over  the  entire  acoustic 
domain  of  interest. 

Figure  5  shows  the  amplitude  response  of  a  uniaxial  velocity  sensor  aligned  along  the 
y-axis,  or  the  amplitude  Vy  of  the  velocity  component  Vy. 
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A  cross-section  of  this  element  response,  showing  the  cosine  angular  dependence,  is 
given  in  figure  6.  The  amplitude  of  Fy  is  reduced  by  one-half  at  ±60°.  Responses  for  Vx  and 
are  shown  in  figures  7  and  8.  For  illustration,  both  of  the  amplitude  lobes  of  Vx  and  V^  are 
shown  positive;  the  polarity  is  suppressed.  Figures  9, 10,  and  1 1  show  the  element  response 
combinations  possible  from  a  biaxial  vector  sensor.  In  figure  9,  the  vector  magnitude  of  the 
amplitudes  of  Fx  and  Fy  are  shown.  Notice,  in  figure  10,  that  a  biaxial  sensor  which  measures 
only  velocities  in  the  xz-plane  produces  a  null  at  broadside  steerings  (0  =  0°). 

The  element  response  for  a  triaxial  vector  sensor  is  given  in  figure  12,  and  as  with  the 
half-omnidirectional-pressure  sensor,  the  element  response  is  uniform  for  all  elevation  and 
azimuth  arrival  angles.  Hence,  the  magnitude  of  a  triaxial  vector  sensor  yields  an  element 
response  equivalent  to  a  scalar  pressure  sensor.  Some  may  then  question  what  benefit,  in  terms 
of  element  sensitivity,  triaxial  vector  sensors  have  over  pressure-sensors.  One  notable  benefit  is 
described  in  the  following  section,  namely,  the  ability  to  steer  a  single  vector  sensor  in  the 
planewave  arrival  direction. 


Figure  6.  Cross  Section  of  the  Amplitude  Response  of  the  vyComponent 
of  a  Uniaxial  Vector  Sensor  versus  Azimuth  Angle  (p-0°) 
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(|),  ELEVATION  ANGLE  .  ELEVATION  ANGLE 


0,  AZIMUTH  ANGLE 


0,  AZIMUTH  ANGLE 


Figure  7.  Amplitude  Response  of  the 
Vx-Component  of  a  Uniaxial  Vector  Sensor 
on  an  Ideal  Pressure  Release  Boundary 


Figure  8.  Amplitude  Response  of  the 
V ^Component  of  a  Uniaxial  Vector  Sensor 
on  an  Ideal  Pressure  Release  Boundary 


0,  AZIMUTH  ANGLE 

Figure  9.  Magnitude  of  the  Vx-  and  ly- 
Components  of  a  Biaxial  Vector  Sensor  on 
an  Ideal  Pressure  Release  Boundary 


0,  AZIMUTH  ANGLE 

Figure  10.  Magnitude  of  the  Vx-  and 
Components  of  a  Biaxial  Vector  Sensor  on 
an  Ideal  Pressure  Release  Boundary 
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0,  AZIMUTH  ANGLE 

Figure  11.  Magnitude  of  the  Vy-  and 
Components  ^  a  Biaxial  Vector  Sensor 
on  an  Ideal  Pressure  Release  Boundary 
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Figure  12.  Magnitude  of  the  v^-,  Vy-,  and 
v^Components  of  a  Triaxial  Vector  Sensor 
on  an  Ideal  Pressure  Release  Boundary 


BEAMFORMING  A  SINGLE  VECTOR  SENSOR 

It  is  possible  to  form  a  broad  acoustic  beam,  which  is  steered  in  the  direction  of  the 
incident  wavevector,  k ,  using  a  single  triaxial  vector  sensor.  A  biaxial  sensor  may  be  steered  as 
well,  but  only  in  one  angle,  such  as  azimuth. 

Consider  the  vector-sensor  response  given  in  equation  (18),  with  the  sensor  located  at  the 

origin. 


v^^^(e,(|),0)  =  2{cos(0^)n;^  +cos(0^)n^  +cos(02)n2}.  (19) 

Then,  as  outlined  in  reference  4,  the  dot  product  of  the  normalized  velocity  field  vector  at  the 
angle  at  which  the  element  is  steered  (65,  <j)j),  and  the  velocity  vector  at  arrival  angle  (0,  (j)),  may 
be  taken,  yielding 


v<^>(e,,(|),,o)-v(^>(e,(|),o) 


v(3). 


v'-(e„(t)„o) 


v(3>(e,<|),o) 


=  /^\e„(t)„0,(|)) 

=  cos((l)^)cos((|))cos(0-05)  +sin((t)j)sin((t)). 


(20) 


The  product  g^^^(05,  (t)j,  0,  <}))  reaches  a  maximum  of  unity  when  (0^,  (l)^)  =  (0,(1)).  For  a  biaxial 
sensor,  measuring  the  components  Vx  and  Vy,  the  steereci  element  response  is 


90 
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<().  ELEVATION  ANGLE 


g'^>(e,.9)  =  cos(e-6.). 


(21) 


Figures  13  and  14  present  the  amplitude  response  and  the  associated  beam  pattern  (log  of 
the  squared  amplitude  response)  of  a  single  triaxial  vector  sensor  steered  to  45°  azimuth  and  30° 
elevation.  The  amplitude  is  not  symmetrical  about  the  steered  direction,  i.e.,  the  beamwidth  is 
broader  in  azimuth  than  elevation.  This  is  caused  by  the  cosine  directionality  of  each  component 
amplitude,  the  summed  response  being  a  combination  of  cosines. 

Shown  in  figure  15  are  beam  patterns  for  a  triaxial  element  steered  to  0°  azimuth  and 
various  elevation  angles.  In  elevation,  the  power  in  each  element  beam  is  reduced  by  one-half  at 
+60°,  as  predicted  for  cosine  directivity  and  shown  previously.  The  half-power  beamwidth  does, 
however,  vary  with  other  element  steerings.  For  example,  figure  14  illustrates  the  beam  pattern 
for  the  amplitude  response  given  in  figure  13.  Here,  the  half-power  beamwidth  occurs  when  the 
arrival  angles  are  such  that  g(3^(45°,  30°,  0,  (|))  equals  one-half.  All  of  the  element  beamwidths  are 
frequency-independent. 


e,  AZIMUTH  Angle  e,  azimuth  angle 


Figure  13.  Amplitude  Pattern  of  a  Single  Figure  14.  Beam  Pattern  of  a  Single 

Triaxial  Vector  Sensor  Steered  to  Triaxial  Vector  Sensor  Steered  to 

45°  Azimuth  and  30°  Elevation  45°  Azimuth  and  30°  Elevation 


With  additional  sensors,  the  accuracy  of  the  estimation  of  the  signal  arrival  angle  from 
steering  individual  sensors  would  improve  -  assuming  the  steering  errors  at  each  sensor  were 
uncorrelated.  Hence,  one  could  consider  steering  a  horizontal  line  array  in  elevation  by  averaging 
over  each  single-element  beam  response.  The  averaging  process  would  provide  little  resolution 
improvement,  so  the  directivity  in  elevation  would  be  limited,  as  indicated  in  figure  15.  The 
elements  in  the  horizontal  direction  can,  of  course,  be  phase  delayed,  weighted,  and  summed, 
which  would  provide  azimuthal  resolution.  Thus,  one  possible  array  beamforming  procedure  is 
to  steer  each  vector-sensor  element  in  the  direction  of  Ihe  estimated  arrival  angle  and  beamform 
the  resulting  responses  conventionally. 
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Figure  15.  Beam  Pattern  in  Elevation  of  a  Single  Triaxial  Vector  Sensor 
Steered  to  Elevation  Angles  of  0°,  45°,  and  90° 


BEAMFORMING  AN  ARRAY  OF  VECTOR  SENSORS 

The  line  array  depicted  in  figure  1  consists  of  ten  elements,  with  interelement  spacing  dx- 
This  array  configuration  allows  for  element-to-element  phase  variations  along  the  x-axis  only. 
That  is,  for  a  fixed  incident  planewave,  the  phase  at  each  element  is  constant  with  respect  to  the 
y-  and  z-axis  and  varies  only  with  respect  to  the  x-axis.  Therefore,  phase  can  be  delayed  (or  the 
array  can  be  steered)  only  by  delaying  phase  components  along  the  x-axis,  such  as  is  the  case 
with  scalar  pressure  sensors.  This  is  fhndamental;  with  conventional  array  beamforming, 
regardless  of  the  type  of  array  (line,  planar,  or  volumetric),  the  phase  delays  are  no  different  for  a 
scalar  or  vector-sensing  array. 

A  possible  approach  to  beamforming  an  array  of  vector  sensors  would  be  to  delay, 
weight,  and  sum  the  velocity  components  from  each  triaxial  vector  sensor  separately.  Thus,  for 
the  ideal  baffle  configuration 

iv  _  _  _ 

(|))  =  2  sin(e)  cos((l))^  +2  cos(0)  cos((|))2]  fiy 

«=1  n=\ 


N 

«=i 


(22) 
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where  the  amplitude  shading  coefficients,  are  arbitrary.  The  wavevector  k , 

defined  by  equations  (1 1),  (12),  and  (13),  corresponds  to  the  arriving  acoustic  planewave,  and 
is  the  wavevector  to  which  the  array  is  steered.  It  is  also  noted  that  the  exponential  functions  in 
equation  (22)  may  be  written  more  explicitly  as 


^{k-k,yr„ 


(0 


i  7  i  <1>)  +  <|))  +  ^ 


„r(e.(i.)} 


(23) 


where 

a(0,(j))  =  cos((t>)sin(9)  -  cos((t)^)sin(0j),  (24) 


p(0,(t))  =  cos((|))cos(0)  -  cos((t>^)cos(0^),  (25) 

and 

y(0,(|))  =  sin((|))-sin((|),).  (26) 


One  may  then  be  tempted  to  simply  take  the  norm  (or  length  of  the  weighted  sum  of  all 
the  vector  components)  of  in  equation  (22)  as  a  means  to  process  the  vector  components. 
That  is,  the  norm  may  be  written  as 


N 


II I  ^xn 

y  V«=l 


2  r 


N 


V«=l 


yn 


r 

+ 


J 


N 
V^=l 


zn 


where  exp{/(^  -k)-r^} ,  etc.  Clearly,  this  approach  results  in  taking  the  square-root  of 

the  sum  of  the  squared  vector  velocity  components;  this  is  a  nonlinear  ad  hoc  processing 
technique  that  should  be  avoided. 


Another  approach  to  array  beamforming,  as  previously  described,  would  be  to  steer  each 
vector  sensor  in  equation  (22)  towards  the  source,  and  then  obtain  the  squared  magnitude  of  the 
resulting  expression.  That  is,  form  the  inner  product  of  the  summed  velocities,  ,  with  a 
steered  vector  of  unit  amplitude  and  take  the  magnitude-squared  of  the  product.  If  it  is  assumed 
that  the  weighting  of  each  component  of  every  vector  sensor  is  identical,  i.e.. 


w, 


=  wi,  '  =  w, 


n  5 


and,  for  the  line  array  in  figure  1,  =  [x„,0.,0.],  then  this  approach  yields 
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s'”(e,(t.)=||v?(e,(t.)-v«>(e„^>.)f 


Z'*.«p(i-j:,a(e,it)))| 


Similarly, 


BP)(e,.|,)=4{g»>(e„e)cos«)}2 


«=1 


(27) 


(28) 


There  is  no  equivalent  expression  for  a  uniaxis  vector-sensing  array  since  such  elements 
cannot  be  steered  in  azimuth  or  elevation.  Therefore,  the  uniaxial  beam  response  (for  measuring 
the  vector  component  v^,)  would  be  written  simply  as 


v^^\e,(t)) 

2  0  0 

N  ^ 

=  4cos^(0)cos^((j)) 

Xw'„exp(/-x„a(0,(f))) 

n-\ 

(29) 


As  an  aside,  it  is  noted  that  the  summations  given  above  can  be  written  in  a  more 
compact  form  for  certain  array  configurations.  For  example,  with  uniform-shading  weights,  and 
equal  element  array  spacing,  Lagrange's  trigonometric  identity  gives 

N  ,  . 

2]  yvn  exp^zy  x„a(0,  (t))j  =  sin(iV/i )  /  sin(44)  if  xj  =  -X]t«j ,  (30) 

n=\ 

where 

A  =  7t-^a(0,<|)). 


This  second  approach  to  beamforming,  however,  may  be  limited  in  its  capability.  The 
restriction  of  weighting  values  to  be  essentially  the  direction  cosines  and  constant  from  element 
to  element  disallows  emphasis  of  some  element  outputs  over  others  and  it  imposes  the  same 
relative  weighting  to  the  three  velocity  components. 

A  more  general  linear  processing  technique  would  form  the  inner  product  of  the  velocity 
at  each  element  location  with  an  arbitrary  element  weighting  vector,  ,  which  now  may  vary 
at  each  element  array  location.  The  resulting  scalar  field  may  then  be  delayed  and  summed  in  a 
conventional  manner.  Thus,  in  general, 


n-\ 


(31) 
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|2 


«=l 


(31) 


and  for  the  line  array, 


|2 


+  exp(/^x„  a(0,<^)) 


(32) 


Hence,  the  result  here  is  a  linearly-weighted  combination  of  all  three  measured  velocity 
components.  That  is,  this  approach  yields  the  square  of  the  sum  of  the  weighted  velocity 
components.  A  special  case  of  the  general  weighting  vector,  ,  would  be  constant  direction 

cosine  weights. 


The  most  general  linear  processing,  likely  only  realizable  in  a  free-field  environment, 
would  be  to  augment  equation  (M)  with  a  measurement  of  acoustic  pressure.  Processing  all  4N 
measured  quantities,  i.e.,  pressure  and  three  velocity  components  at  each  element  location,  yields 


N 


n=l 


(33) 


There  are  now  4N  degrees  of  freedom  to  manipulate,  in  order  to  maximize  some  measure 
of  performance,  such  as  the  array  ouq)ut  signal-to-noise  ratio.  These  4N  degrees  of  freedom  are 
the  N  pressure  weights  {w^}  and  the  N  velocity  weights  for  each  orthogonal  direction,  namely, 

{  }>  { }  for  1  <  n  <  N.  This  increased  number  of  degrees  of  freedom  should 

allow  for  considerably  improved  performance  relative  to  pressure-alone  processing. 

For  reference,  the  power  response  for  a  line  array  of  baffled  scalar  pressure  sensors  is 


#^\0,(t))  =  4 


N 


Xi^«exp(ijx„a(0,<t)))| 
«=1 


(34) 


Hence,  from  equation  (27),  it  is  seen  that  the  acoustic  power  from  a  uniformly- weighted 
line  array  of  pressure  sensors  is  equivalent  to  that  from  an  array  of  triaxial  vector  sensors  at  the 
array's  maximum  response  axis.  At  elevation  angles,  <|)  =  ±kI2,  the  uniaxial  (equation  (29))  and 
biaxial  (equation  (28))  line  arrays  considered  here  receive  no  power,  and  the  uniaxial  vector¬ 
sensing  array  also  has  a  null  response  at  azimuthal  endfire.  Such  is  not  the  case  for  the  half- 
omnidirectional  pressure-sensing  array. 
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NUMERICAL  RESULTS 


COMPARISON  OF  BEAM  PATTERNS 

The  beam  pattern  of  an  ideally-baffled  10-element  triaxial  vector-sensing  line  array,  with 
half- wavelength  element  separations  and  steered  to  broadside,  is  shown  in  figure  16,  which 
employs  a  gray  scale  from  -30  dB  to  0  dB.  For  comparison,  figure  17  illustrates  the  beam  pattern 
from  a  pressure-sensing  array  steered  to  broadside.  For  broadside  steering,  the  beam  patterns  from 
both  array  types  are  alike.  The  triaxial  vector  sensors  provide  greater  vertical  directivity,  due  to 
the  cosine  dependence,  than  the  half-omnidirectional  pressure  sensors.  Figure  18  compares  the 
primary  sidelobe  structure,  in  a2dmuth,  of  both  arrays.  Again,  the  cosine  angular  dependence  of 
the  vector  sensor  lowers  the  sidelobe  levels,  the  reduction  increasing  as  the  azithumal  arrival 
angle  increases.  This  overall  sidelobe  reduction  produces  an  insignificant  broadening  of  the 
azimuthal  mainbeam  width. 


0,  AZIMUTH  ANGLE  6,  AZIMUTH  ANGLE 


Figure  16.  Beam  Pattern  at  Broadside  of  a  Linear  Figure  1 7.  Beam  Pattern  at  Broadside  of  a  Linear 
10-Element  Triaxial  Vector-Sensing  Array  10-Element  Pressure-Sensing  Array 

Having  Uniform  Shading  Weights  Having  Uniform  Shading  Weights 


The  differences  in  vector-  and  scalar-array  patterns  become  more  apparent  when  the  array 
is  steered  to  oblique  angles  of  incidence.  For  example,  figure  19  illustrates  the  beam  pattern 
produced  from  a  triaxial  vector-sensing  array  steered  to  45°  azimuth  and  30°  elevation. 

Referring  back  to  figure  14,  it  is  clear  that  the  array’s  beam  pattern  is  the  product  of  the  steered 
element  response  and  the  array's  response.  The  beam  pattern  from  a  pressure-sensing  array  is 
presented  in  figure  20.  In  both  figures,  the  null  structure  of  the  array's  response  is  identical,  as 
indicated  by  equations  (27)  and  (34).  The  differences  seen  are  due  to  the  inherent  vector-sensor 
directionality  and  single-element  steering  of  each  vector  sensor.  Therefore,  at  oblique  angles  of 
incidence,  a  vector  sensor  may  provide  an  increase  in  array  directivity  over  that  of  a  pressure- 
sensing  array.  This  increase  is  quantified  in  the  next  section. 
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(]),  ELEVATION  ANGLE 


Figure  19.  Beam  Pattern  of  a  Linear 
10-Element  Triaxial  Vector-Sensing  Array 
Having  Uniform  Shading  Weights  Steered  to 
45°  Azimuth  and  30°  Elevation 


Figure  20.  Beam  Pattern  of  a  Linear 
1 0-Element  Pressure-Sensing  Array  Having 
Uniform  Shading  Weights  Steered  to 
45°  Azimuth  and  30°  Elevation 


17 


FREE-FIELD  DIRECTIVITY  CALCULATIONS 


A  comparison  of  free-field  array  directivity  of  a  line  array  of  scalar  sensors  detecting 
acoustic  pressure  and  an  array  of  vector  sensors  detecting  acoustic  particle  motion  in  up  to  three 
orthogonal  directions  is  presented  here.  Appendix  A  includes  an  analytical  derivation  for  the 
directivity  of  some  vector-sensing  line  arrays.  The  calculations  presented  here,  though,  were 
based  on  the  derivative  given  in  Appendix  B.  This  derivation  maximizes  directivity  by  choosing 
array  shading  weights  constrained  to  be  real-valued,  which  maximize  the  array  signal-to-noise 
ratio  assuming  a  surrounding  spherically  isotropic  noise  field. 

Figure  21  compares  the  maximum  directivity  index  of  an  N-element  pressure  sensing  line 
array  to  an  uniaxial  vector-sensing  array  which  measures  only  the  y-component  of  acoustic 
particle  velocity.  The  elements  of  the  arrays  are  equispaced  at  half-wavelength  spacing  and  are 
aligned  along  the  x-axis  in  the  fi:ee  field.  The  arrays  have  been  steered  to  broadside  (G^  =  0°) 

in  figure  21  and  endfire  (Gj  =  90°,  (jij  =  0°)  in  figure  22.  Also  shown  in  the  figures  is  the 
directivity  index  for  an  array  which  measures  both  acoustic  particle  velocity  and  pressure. 

For  a  planewave  signal  arriving  at  broadside,  that  is,  along  the  y-axis,  there  would  be  no 
particle  velocity  in  the  x-  or  z-directions.  Hence,  for  this  array  steering,  there  would  be  no 
additional  directivity  gains  from  measuring  the  velocities  Vx  or  v^.  Similarly,  for  planewaves 
arriving  at  endfire,  only  the  x-component  of  acoustic  particle  velocity  is  non-zero. 

Since  the  directivity  calculations  are  based  on  a  half- wavelength  equispaced  line,  a  10- 
element  pressure-sensing  array  has  a  length  of  5X,  for  a  directivity  of  10  dB.  Notice,  in  figure  21, 
that  the  imiaxial  vector  sensing  array  achieves  10  dB  of  directivity  with  half  that  number  of 
elements  and  has  a  length  of  2.5X.  The  array  that  measures  both  pressure  and  velocity  provides 
10  dB  of  directivity  at  a  length  of  only  1.5X,  and  the  total  number  of  channels  (i.e.,  6)  is  still  less 
than  the  niunber  of  elements  on  the  pressure-sensing  array. 

The  differences  between  the  directivity  gains  shown  in  these  figures  is  due  to  the  orienta¬ 
tion  of  the  vector  sensors  with  respect  to  the  incident  planewave.  In  figure  21,  the  difference  in 
directivity  between  pressure  sensing  and  vector  sensing  (v_y-component  only)  remains  essentially 
constant  with  increasing  array  aperture.  This  is  because  the  azimuthal  beamwidth  narrows,  with 
increasing  aperture,  identically  for  the  pressure  and  vector  sensors.  Hence,  the  gain  shown  is  due 
to  a  constant  difference  between  elevation  beamwidths.  The  uniaxial  vector  sensor  has  a  cosine 
element  response,  whereas  the  pressure  sensor  is  assumed  omnidirectional.  In  elevation,  the 
directivity  of  the  vector  sensor  is  twice  that  of  the  pressure  sensor.  Hence,  there  is  a  constant 
difference  of  3  dB. 

In  figure  22,  the  difference  between  the  directivity  index  of  the  vector-sensing  array  and  a 
pressure-sensing  array  diminishes  with  increasing  array  aperture.  Again,  the  endfire  beamwidth 
narrows  identically  for  both  the  pressure-  and  vector-sensing  arrays  with  increasing  aperture. 
However,  since  the  arrays  are  steered  to  endfire,  the  effect  of  the  broad  angular  difference 
between  the  cosine  and  omnidirectional  element  response  becomes  negligible  in  comparison  to 
the  endfire  beamwidth.  That  is,  with  increasing  aperture,  the  directivity  is  determined  entirely  by 
the  array  beamwidth  and  is  minimally  influenced  by  the  individual  element  response.  This  is  a 
feature  unique  to  endfire  steerings. 
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DIRECTIVITY  INDEX  [dB]  DIRECTIVITY  INDEX  [dB] 


1 .5  X  2.5  ^  5  ^  [LENGTH] 
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N  [NUMBER  OF  ELEMENTS] 


Figure  22.  Maximum  Directivity  for  an  N-Element  Equispaced  Line  Array  at 
Endfire  (6$  =  90°)  Assuming  Optimal  Real  Weights 
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Also  shown  is  the  directivity  of  an  array  that  measures  both  acoustic  pressure  and  particle 
velocity.  Including  the  pressure  increases  the  directivity  for  both  array  steerings  and  the 
increased  gain  appears  to  be  directly  proportional  to  the  number  of  array  elements,  N.  That  is,  the 
directivity  factor  can  be  expressed  as  DFp^  =  DF^  +  N,  for  optimal  weightings  at  these  steerings. 
This  simple  relationship  does  not  appear  to  hold  for  general  array  steerings,  as  will  be 
demonstrated  in  figures  23  through  25.  Notice  that,  for  a  single  element  (N  =  1)  and  for  either 
array  steerings,  the  uniaxis  vector  sensor  provides  a  4.8  dB  gain  over  the  omnidirectional 
pressure-sensing  element.  Measurement  of  both  pressure  and  velocity  yields  a  6  dB  gain. 
Appendix  A  provides  a  description  of  the  source  for  these  gains. 

In  figure  23,  a  comparison  is  made  between  the  variation  of  directivity  with  azimuth  for  a 
pressure-sensing  and  vector-sensing  line  array.  A  10-element  line  array,  again  aligned  along  the 
x-axis  and  in  the  free  field,  is  assumed  and  the  elevation  angle  is  fixed  at  broadside.  In  figure  24, 
the  elevation  angle  is  fixed  at  15°,  whereas  in  figure  25,  the  elevation  angle  is  30°.  In  the  free 
field,  the  optimal  directivity  of  an  array  of  omnidirectional  pressure  sensors  is  proportional  to  the 
number  of  array  elements,  N,  for  all  array  steerings.  Hence,  the  constant  directivity  index  of  10 
dB  for  the  pressure-sensing  array  is  observed.  Notice  in  figure  23  that  the  array  of  uniaxial  vector 
sensors,  denoted  Vy,  provides  3  ^  gain  over  pressure  sensors  for  most  azimuthal  steering  angles. 
This  gain  is  due  to  the  differences  in  the  cosine  and  omnidirectional  element  response  patterns. 
However,  at  oblique  elevation  angles,  this  gain  appears  to  diminish  in  proportion  to  the  cosine  of 
the  elevation  angle. 

Biaxial  and  triaxial  sensors  further  improve  directivity,  particularly  for  incident  angles  far 
from  broadside  where  the  gains  can  be  very  large.  For  example,  in  figure  24  the  directivity  index 
of  the  uniaxial  array  is  5  dB  at  77°  azimuth.  On  the  other  hand,  the  triaxial  array  directivity  is 
14  dB,  providing  a  directivity  gain  of  9  dB  against  isotropic  noise.  Notice  that  the  maximum 
gain  is  achieved  when  the  pressure  and  all  three  components  of  particle  velocity  are  measured 
and  that  this  optimal  gain  is  constant  for  all  array  steerings. 


Figure  23.  Directivity  Versus  Azimuthal  Angle  (Os)  for  a  10-Element 
Equispaced  Line  Array  at  Elevation  Angle  (4>s~^°) 
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Figure  24.  Directivity  Versus  Azimuthal  Angle  (6^  for  a  10-Element 
Equispaced  Line  Array  at  Elevation  Angle  =  15°) 


Figure  25.  Directivity  Versus  Azimuthal  Angle  (dg)  for  a  10-Element 
Equispaced  Line  Array  at  Elevation  Angle  (^g  =  30°) 
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THE  EFFECT  OF  ARRAY  CURVATURE 


Figure  26  depicts  a  section  of  a  vertical  line  array,  curved  about  the  x-axis,  with  elements 
lying  in  the  yz-plane.  For  this  configuration,  the  component  normals  ity'  and  vary  with 
element  position  and  may  be  written  as 

fiy.  =  cos(p„)n^  +  sin(P„)«^  (35) 

and 

=  -  sin(p„)«^  +  cos(P„)n2,  (36) 

where  p„  is  the  angle  formed  between  the  n-th  element  normal  and  the  y-axis. 

Now,  the  cosine  projection  angles  for  the  y-  and  z-velocity  components  (as  in  equation 
(14),  taking  the  dot  product  of  the  signal-arrival  vector,  k ,  with  the  component  normals)  become 

cos(©y)  =  cos((j))cos(e)cos(p„)  +  sin(<|))  sin(p„),  (37) 

and 

cos(02<)  =  sin((j))cos(P„)  -  cos((t))  cos(0)sin(P„).  (38) 


k 


\ 


1 

t 

Figure  26.  Section  of  a  Curved  Vertical  Line  Array 
Comprised  of  Biaxial  Vector  Sensors  (vy,  v^)  and 
Component  Normals  (fiy’ ,  n^’) 

Using  the  above  result,  the  signal  gain  of  a  curved  line  array  of  uniaxial  vector  sensors, 
which  measures  the  vector  component,  Vy,  to  an  array  of  biaxial  sensors,  measuring  both  the  v^;- 
and  Vz-components,  may  be  compared.  These  gains,  which  are  independent  of  frequency,  reveal 
how  curvature  and  vector  element  responses  combine  and  affect  array  performance.  The  uniaxial 
and  biaxial  baffled  array  sensitivities  can  be  defined,  for  all  steerings,  respectively,  as 
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^  cos(©y  (0, ,  (|)  jV„(<t)  J 


n=l 


(39) 


M^(e„<t.,)  =  4 


N 


n=l 


2]  cos2(0^,(e„(t),))+cos2(©,,(e,,(j,,))U  ((^,) 


(40) 


where  it  is  recognized  that  the  n-th  element  shading  weight,  >v„(<|)5),  may  be  a  function  of  position 
and  elevation  steering  angle.  The  biaxial  array  sensitivity,  Myz,  is  arbitrarily  assumed  to  have 
equivalent  shading  of  the  and  velocity  components. 

For  a  given  signal  power  level,  the  signal  gain  resilized  from  a  curved  array  of  biaxial 
vector  sensors,  over  that  of  uniaxial  sensors,  is  defined  as 


5G(e„(t>,)  =10Zngio 


^yi^sAs) 


(41) 


For  broadside  azimuth  steering,  the  expression  for  signal  gain  simplifies  to 
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|2 


l^«(<t>5) 


5G(0,(t)5)  =10Zogloi 


\n=l  I 
N 

Xcos(<|)^-p„)w„((|)J 


\n=l 


(42) 


The  signal  gain  of  a  fully-populated  cylindrically-curved  line  array,  steered  to  broadside, 
is  given  in  figure  27.  The  biaxial  sensors  improve  the  conformal  array's  signal  gain  at  all 
elevation  steerings,  providing  a  minimum  gain  of  approximately  3  dB  at  ±45'  and  a  maximum 
gain  of  nearly  4  dB  at  both  broadside  and  endfire  elevation  angles. 
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Figure  27.  Signal  Gain  of  a  Fully-Populated  Cylindrically  Curved  Vertical  Line  Array 
of  Biaxial  Vector  Sensors  with  Uniform  Shading  Weights 


CONCLUSIONS 

Presented  in  this  report  are  formulations  for  the  beam  response  and  directivity  of  arrays  of 
uniaxial,  biaxial,  and  triaxial  vector  sensors.  To  simplify  the  development,  each  orthogonal 
component  of  a  multiaxis  vector  sensor  was  assumed  to  be  collocated  and  the  orientation  of  the 
components  was  unvarying  from  sensor  to  sensor.  However,  the  amplitude  shading  coefficients 
for  each  vector  component  of  particle  velocity  was  allowed  to  be  arbitrary.  The  report  has 
investigated  the  response  of  the  sensors  in  both  the  free  field  and  on  an  ideal  baffle. 

Of  the  results  presented  in  this  report,  the  most  significant  are  the  increased  gains  in 
directivity  an  array  of  vector  sensors  provides,  shown  in  figures  21  through  25,  over  that  of  a 
pressure-sensing  array  in  the  free  field  Directivity  calculations  for  an  ideally-liaffled  line  array 
is  an  on-going  effort  and  have  not  been  presented  here.  Figure  23  indicates  that  a  10-element 
line  array  of  triaxial  vector  sensors,  with  elements  separated  by  one-half  an  acoustic  wavelength, 
has  a  directivity  index  5  dB  greater  than  that  of  an  identical  line  array  of  pressure  sensors  for  all 
azimuthal  array  steerings  between  12°  to  67°.  Hence,  to  achieve  equivalent  array  gain,  the 
acoustic  aperture  of  the  pressure-sensing  line  array  must  be  almost  four  times  greater  than  that  of 
a  vector-sensing  array.  This  implies  that  vector-sensing  arrays,  albeit  with  a  larger  niunber  of 
output  channels,  may  be  physically  smaller  than  pressure-sensing  arrays,  Avithout  a  loss  in 
directivity.  With  pressure  sensors,  the  only  way  to  increase  directivity  is  to  increase  array 
aperture.  Therefore,  for  pressure-sensors,  gains  in  directivity  are  fundamentally  limited  by  the 
size  of  the  array.  However,  for  a  fixed  array  aperture  and  a  tlu-ee-fold  increase  in  chaimel  count, 
significant  directivity  gains  may  be  realized  with  an  array  of  triaxial  vector  sensors.  Increasing 
the  number  of  processing  channels  does  add  complexity  and  expense  to  an  array  design  but  such 
design  challenges  are  modest  compared  to  trying  to  increase  array  aperture. 
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Another  interesting  result  described  here  is  the  ability  to  form  a  broad  acoustic  beam, 
steered  in  the  direction  of  an  incident  signal,  using  a  single  multiaxis  vector  sensor.  Hence,  a 
vector-sensing  line  array,  oriented  say,  along  the  x-axis  as  shown  in  figure  1,  will  provide 
resolution  of  both  the  azimuth  and  elevation  arrival  angles.  It  was  noted,  however,  that  the 
resolution  of  the  elevation  angle  would  be  limited  to  a  beamwidth  proportional  to  the  cosine  of 
the  elevation  angle.  That  is,  &e  beamwidth  is  ±60°  at  the  point  at  which  the  power  relative  to  the 
maximum  response  axis  is  down  by  one-half 

Figures  4  and  5  illustrated  the  angular  sensitivity  of,  respectively,  an  individual  pressure 
sensor  and  the  y-component  of  a  uniaxial  vector  sensor.  By  definition,  the  uniaxi^  vector  sensor 
element  response  is  more  directional,  and  it  is  this  directionality  which  provides  the  constant  3 
dB  increase  in  the  array  directivity  index  shown  in  figure  21 .  Specifically,  the  cosine  element 
response  has  an  element  directivity  that  is  twice  that  of  an  omnidirectional  pressure  sensor.  This 
gain  then,  for  broadside  array  steering  only,  is  due  simply  to  the  assumed  differences  in  the 
pressure  and  uniaxial  element  responses. 

It  must  be  noted  that  baffled  pressure  sensors  often  have  an  element  response  pattern 
approximated  well  by  a  cosine  raised  to  a  fractional  power  approaching  unity.  For  these 
configurations,  the  directivity  index  for  a  linear  pressure-sensing  array  and  any  type  of  linear 
vector-sensing  array  steered  to  broadside,  would  be  approximately  equivalent.  A  vector-sensing 
array  will  provide  a  constant  directivity  gain  of  3  dB  over  a  pressure-sensing  array  when  it  is 
assumed  the  pressure  sensors  are  half  or  omnidirectional.  However,  additional  gains  are  seen,  as 
shown  in  figures  23, 24,  and  25,  for  other  array  steerings.  Here,  the  pressme-sensing  directivity 
remains  constant,  whereas  the  multiaxis  vector  sensor  gains  vary  wiA  array  steering  and  these 
gains  are  due  to  measuring  the  additional  vector  components. 

A  brief  investigation  of  curvature  issues  indicates  that  multiaxis  vector  sensors  reduce 
signal  Josses  due  to  highly-curved  array  apertures.  An  example  is  given  of  a  fully-populated 
curved  vertical  line  array  comprised  of  eitihier  biaxial  vector  sensors  or  uniaxial  vector  sensors. 
The  biaxial  vector  sensors  improved  the  signal  gain,  independent  of  frequency,  over  a  uniaxial 
vector  sensing  array  by  as  much  as  4  dB.  Similar  gains  would  be  seen  in  a  comparison  to  a 
pressure-sensing  array  and,  in  general,  multiaxis  vector  sensors  will  improve  a  conformal  array's 
signal  gain  for  all  array  steerings. 

Finally,  it  has  been  shown,  particularly  in  figures  24  and  25,  that  it  is  desirable  to  have  an 
array  which  measures  pressure  in  addition  to  the  three  components  of  acoustic  particle  velocity. 
Optimal  processing  of  arrays  which  measure  all  four  quantities,  pressure  and  the  three 
components  of  velocity,  is  the  subject  of  future  research. 
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APPENDIX  A 


ANALYTICAL  EVALUATION  OF  DIRECTIVITY 
FOR  HALF-WAVELENGTH  ARRAY  SPACING 

Array  directivity  is  defined  as  the  ratio  of  the  output  signal-to-noise  ratio  (SNR)  of  an 
array  to  the  input  SNR  at  an  omnidirectional  element  in  a  spherically  isotropic  noise  field. 
Directivity  is  a  fundamental  measure  of  array  performance;  the  directivity  factor,  DF,  of  an  array 
in  the  firee  field  can  be  expressed  as^ 

47iR(e„f,)  (A-1) 

~  Inn 

J  jR(e,<|))sin((t))4J0 

0  0 

where  R(0,<t))  is  the  array’s  beam  pattern.  The  directivity  index  is  10  log(DF).  Note  that  to 
simplify  the  following  derivations,  the  elevation  angle  ^  will  now  be  measured  from  the  vertical 
z-axis  and  the  azimuthal  angle  0  will  be  measured  from  the  x-axis  in  the  direction  of  the  y-axis. 
Also,  the  line  array  elements  are  now  equispaced  along  the  z-axis  at  spacing  d=  X/2. 

Calculation  of  the  DF  typically  requires  the  numerical  evaluation  of  the  above  double 
integral.  However,  for  certain  array  geometries  and  fi’equencies  (e.g.,  a  frequency  at  which  an 
acoustic  wavelength  equals  tvdce  the  interelement  line  array  spacing,  =  XU),  the  integrals  may 
be  evaluated  exactly.  These  frequencies  correspond  to  the  array's  upper  aliasing  frequency  limit 
and  maximum  array  directivity. 

The  maximum  array  response  at  steering  angles  (0s,(t)s)  can  be  obtained  from  equation 
(33).  In  the  free  field  and  with  elements  placed  along  the  z-axis  at  one-half  an  acoustic 
wavelength,  the  maximum  array  responses  may  be  written  as 
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and,  for  triaxial  vector  sensors, 


(A-5) 


N-l 


«=0 


(A-6) 


where  a(0s,<t)s)  =  cos(0s)  sin((|)s),  6(0s,<t>s)  =  sin(0s)  sin((|)s),  and  c((j)s)  =  cos((t)s). 

The  effort  in  obtaining  expressions  for  directivity  factors  then  is  in  analytically  evaluating 
the  double  integral  given  in  equation  (A-1). 


DIRECTIVITY  FACTORS  FOR  A  SINGLE  ELEMENT,  N  =  1. 

It  is  straightforward  to  derive  the  directivity  factor  of  a  single  omnidirectional  pressure 
sensor.  Namely, 


2n  n 
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exp|mJt(cos(<t))  -  cos((()^  )| 
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sin((|))  dQdi^  =  4it(  .  (A-7) 


Hence,  upon  substituting  equation  (A-7)  and  (A-2)  into  equation  (A-1), 


477(4'’*)" 


(A-8) 


Therefore,  the  directivity  index  is  zero  and  the  single  pressure  sensor  has  no  directivity. 

A  similar  calculation  for  an  uniaxial  vector  sensor  that  measures  the  V;^  component  of 
particle  velocity,  as  in  equation  (A-3),  gives 


\  J[><*o^^cos(0)sin((t>)J  sin((j))i/0J(t)  =  ^(w^^^)^.  (A-9) 

0  0 

Upon  substitution. 


DF(^x)  =  3cos^(0^)sin^((|)  J. 


(A-10) 


A-2 


(A-11) 


Likewise,  for  the  velocity  components  Vy  and  v^. 


=3sin^(e5)sin^((l)^), 

and 

DF(^)  =3  cos^((t>^).  (A-12) 

Equations  (A- 10),  (A-1 1),  and  (A-12)  indicate  that  a  single  velocity  sensor  may  have  a  maximum 
directivity  factor  three  times  greater  than  an  omnidirectional  pressure  sensor,  or  a  gain  in 
directivity  index  of  4.8  dB.  (Figures  21  and  22  illustrate  this  gain.)  These  results  also  indicate 
that  directivity  factors  of  vector  sensors  are  also  highly  dependent  on  azimuth  and  elevation 
angles  (0s,<t)s). 

The  most  general  type  of  single  array  element  would  be  a  triaxial  vector  sensor 
augmented  with  a  measurement  of  pressure,  as  in  equation  (33).  The  maximum  response  of  such 
an  element  is  given  as 


=  (4'’’  +  +  c('t's)4"0 


(A-13) 


and  the  two-fold  integral  becomes 

2n  71 


0  0 


Hence,  the  directivity  factor  of  a  single  element  that  measures  all  3  components  of  particle 
velocity  and  acoustic  pressure  is 


Fpvi^sAs) 
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(A-15) 


The  goal  now  becomes  determining  the  weights  that  maximize  the  above  ratio.  Without 
loss  in  generality,  one  may  let  =1  and  solve  the  following  equations  simultaneously: 


dDFp^  dDFpy  dDFpy 


(A-16) 


This  leads  to  two  sets  of  optimal  weights,  one  corresponding  to  a  minimum  DF^v  ^nd  the 
other,  given  below,  for  maximum  directivity,  i.e.. 


4^>  =  3*(e„(|.,), 


(A-18) 


4"*=3<;(c|.j).  (A-19) 

Substituting  the  optimal  real  weights  in  equations  (A-17),  (A-18),  and  (A-19)  into 
equation  (A- 15)  yields 

DFp^  =  1  +  3(a(0„(l),)2  +  6(e„<|),)2  +  c((t,,)2j  =  4,  (A-20) 

since  iJf(05,(t)J^  +  =  1  always. 

Therefore,  the  maximum  directivity  for  a  single  sensor  that  measures  all  velocities  as  well 
as  acoustic  pressure  is  DI  =  10  log(4)  =  6  dB,  and  this  holds  for  any  steering. 


DIRECTIVITY  FACTORS  FOR  AN  N-ELEMENT  LINE  ARRAY 

The  integral  given  by  equation  (A-7)  for  pressure  sensors  may  be  integrated  directly  over 
azimuth  angle  0,  and  assuming  real-valued  weights 


N-\  N-l  , 

I  =  2n^  Jexp|/(n  -  w)(cos(<t))  -cos(<|)^))TcJsin((t))d<t).  (A-21) 

^=0  m=0  0 

Letting  u  =  cos((|))  and  integrating  over  u,  results  in  an  exact  expression  for  the  DF  of  a 
line  array  of  pressure  sensors,  i.e., 


~  N-l 


^  N-\ 
V«=0 


n=0 


(A-22) 


This  ratio  is  a  maximum  for  uniform  weights  and  results  in  an  optimal  directivity  of  DI  =  10 
log(N). 


Similar  analytical  expressions  for  the  DF  of  vector-sensing  arrays  are  difficult  to  obtain 
and  require  one  to  approximate  the  DF  by  assuming  a  large  number  of  array  elements  N.  For 
example,  consider  an  array  of  uniaxial  vector  sensors  that  measure  the  Vjj;  component  of  particle 
velocity.  Then 


A-4 


27c  n 

/=  J*  Jsin^(0)sin^(<t)) 
0  0 


\N-\  I 

X  exp[m7i(cos((|))  -  cos((t)^  ))]| 

«=0 


sin(^)4Je.  (A-23) 


Letting  u  =  cos((j))  and  integrating  over  0  yields 

|2 

du. 


1 

N-\ 

1(>- 

u^) 

-1 

w=0 

(A-24) 


Although  this  integral  can  be  evaluated  exactly,  it  is  rather  tedious  and  difficult  to  interpret  and 
approximate.  However,  for  large  N  and  for  Us  -  cos((t)5)  not  close  to  ±1,  one  may  argue  that  the 
exponential  term  in  equation  (A-24)  essentially  samples  (as  in  a  delta  function)  the  integrand  at 
u  =  us.  Therefore, 


-1  «=0  m=0 

(A-25) 

=  27isin^(((>J  j(w^^^)^. 

«=0 


The  DF  then  becomes,  substituting  (A-4)  and  (A-25)  into  (A-1), 


2  sin  ((j)j) 


DFy  = 


V«=0  / 


N-l 

n=0 


(A-26) 


This  agrees  with  the  numerical  calculations  presented  in  figure  21.  This  shows  that,  as  the 
number  of  array  elements  increase  and  for  (j)^  =  90°,  the  directivity  of  a  uniaxial  vector-sensing 
array  approaches  a  gain  that  is  3  dB  greater  than  that  of  a  pressure-sensing  array.  Approximations 
for  DFx,  DFz,  and  even  DFp^  can  be  obtained  in  a  similar  manner  and  will  be  reported  in  a 
subsequent  report. 


A-5/A-6 
Reverse  Blank 


APPENDIX  B 


DERIVATION  OF  MAXIMUM  ARRAY  GAIN 


Consider  an  N-element  array  with  pressure  sensors  and  triaxial  velocity  sensors  at  each 
element  location  Xn,  ^n-  It  Is  presumed  that  each  velocity  sensor  is  oriented  directly  along  its 
corresponding  axis.  If  this  array  is  steered  in  direction  0^,<t)5  by  phase  shifting  the  individual 
element  outputs  Pn(t),  '^xn(0^  '^yn(t)r,  and  Vzn(t),  the  most  general  linearly-processed  array  output 
waveform  is 


N 


«=i 


^pn^pnit)  +  ^xn^xnit)  +  W 


yn 


{t)  +  W 


zn  ^zn 


(B-l) 


where 

Vpn^t)  =  p„{t)l{pc),  k  =  (die, 


e„ (0,  <j))  =  exp|^/^(x„  sin  <j)  cos  0  +  >"„  sin  (|)  sin  0  +  cos  <|))J,  (B-2) 

and  the  4N  weights  {wpn},  {wxn},  {wyn},  {wzn)  are  real,  but  otherwise  arbitrary. 

A  plane-wave  signal  arrival  from  direction  0a,(j)  a  results  in  the  array  output  signal 
waveform 


N 

Y{t)  =  Vi  exp(z©oE4(05.<t>j)e«(ea,'l>a)[wp«  + 

n=l 

(B-3) 

+  sinct)^  cos0a  +  M>yn  sin(j)a  sin0a  +  cos(|)^], 

where  riexp(z©t)  is  the  pressure  waveform  observed  at  the  origin  of  coordinates.  For  notational 
convenience,  define  the  direction  cosines 

=  sin(j)cos0,  Z)^  =  sin(t)sin0,  D^=cos^,  (B-4) 

and  the  two  4Nxl  real  column  matrices 


W=  [Wp]...  WpN  Wxi...  WxN  Wy]...  WyN  Wzi...  WznV, 

D(0,(j))=  [1...  1  Dx...  Dx  Dy...  Dy  Dz...  Dz]T 


(B-5) 


B-l 


Then,  when  the  plane-wave  signal  arrival  angle  coincides  with  the  steering  angle,  namely, 
^  6^5  <l>a  ^  the  array  signal  output  waveform  is  expressible  as 

r,(0  =  F,exp(io)0»'^/)(6„c|.,).  (B-6) 

The  array  signal  output  power  is  then 

Pso  =  |ii('f  =  hf  .  (B-7) 

At  the  same  time,  the  array  output  noise  waveform  can  be  expressed  as 

r„m=wPv^(t),  (B-8) 

the  4Nxl  coluirm  matrix 
'Vpit) 

and  the  four  Nxl  column  matrices 

(B-IO) 

Vy(t)  ^1  (0^j(j)^)  Vjl(f)...^^(0_5  3())^)  V^^(f)|  3 

f^z(0  )  Vzl(0- ••  't>5)  VzAT  (o]^- 

This  leads  to  the  array  output  noise  power  in  the  form 

P«o  =  =  ir^P’r,  (B-ii) 

where  i?  is  a  Hermitian  4Nx4N  covariance  matrix  of  the  totality  of  4N  received  noise  waveforms. 
Because  all  the  weights  in  (B-1)  and  (B-5)  are  real,  (B-1 1)  can  be  written  as 


B-2 


(B-12) 


P^  =  W'^R,W; 

that  is,  the  imaginary  part  of  R  is  not  active  in  Hermitian  form  (B-1 1). 

We  define  the  array  gain  as  the  quotient  of  the  output  signal-to-noise  ratio  to  the  input 
signal-to-noise  ratio.  There  follows 


AG  = 


W^R^W 


(B-13) 


where  we  have  normalized  the  input  noise  pressure  power  at  unity. 

The  4N  optimum  weights  that  maximize  this  array  gain  are 

£>(0^,(1.^),  (B-14) 

while  the  corresponding  optimum  array  gain  is 

=  max^G  =  Z)^(0^,(t)Ji?7'D(0^,<ti5).  (B-15) 

This  relation  was  used  in  the  main  body  of  the  report  to  obtain  the  numerical  values  in  figures  21 
through  25.  The  elements  of  the  real  4Nx4N  matrix  Rf  are  given  in  Appendix  C. 
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APPENDIX  C 


SENSOR  CORRELATIONS  IN  AN  ISOTROPIC  NOISE  FIELD 


In  this  appendix,  ^  is  the  polar  angle  measured  down  from  the  vertical  z-axis,  while  0  is 
the  azimuthal  angle  measured  from  the  x-axis  in  the  direction  of  the  y-axis. 

FUNDAMENTAL  DOUBLE  INTEGRAL 


It  2iz 

Jc/cj)  sin  (j)  I  dQ  exp(za  sin  <t»  cos  0  +  ib  sin  (j)  sin  0  +  ic  cos  (j)) 

0  0 

7t 

=  2ti  ^  sin(t)  exp(zc  cos(|))  sincj)) 

0 

71 

=  271  J  i/(j)  sin(j)  cos(c  cos(|))  Jo((a^  -^b^Y^  sin(|)) 

0 

7C/2 

=  471  j  <i(j)  sin  (j)  cos(c  cos  (j))  ((a^  +  b^Y^  sin  cj)) 

0 

siniQ  ^  /  2  1 2  2V/2 

=  471-^,  g  =  (a  +b^+c)  , 

where  [7;  6.688  2]  is  used  in  the  last  line.  Also  used  is  the  fact  that  sin(c  cos(j))  is  odd  in  (j)  about 
(j)  =  71/2,  whereas  dl  the  other  functions  are  even  in  (j)  about  cj)  =  nil.  This  result  applies  for 
complex  a,b,c  as  well. 

It  is  now  possible  to  evaluate,  in  closed  form,  extensions  of  the  double  integral  above  to 
include  additional  factors  of  the  form 

(sin<j)  COS0)'”  (sin<|)  sin0)"  (cosij))*  , 

where  m,n,k  axe  integers,  by  taking  partial  derivatives  with  respect  to  a  and/or  b  and/or  c.  To  that 
aim,  let 


C-1 


Then, 


^  d  ^  dF  dQ 
F^= — F  = 


cos  Q  sin  Q 


a  =  -a  jiiO)!  Q-, 


°  da  dQ  da 

here  the  term  dQlda  =  alQ  is  used.  The  bracketed  term  is  -1/3  at  Q  =  0.  Also, 

(3-Q^)smQ-3QcosQ 


Fab  =  —  F„  = 


ah  =  ah  jjiQ')!  Q  , 


where  the  bracketed  term  is  1/15  at  0  =  0.  Finally, 


F  =  —  F  = 

wo  a 


[cosg  sing’ 

(3-g^)sing-3gcosg' 

1 — 

(O 

1 

to 

\ _ 

+ 

L  Q'  J 

a 


=  -m)IQ  +  a^  j2{Q)IQ" . 

All  other  partial  derivatives  of  first  and  second  order,  with  respect  to  a  and/or  b  and/or  c,  follow 
from  these  results  by  use  of  symmetry;  that  is,  Q  =  {a^  +  b^  +  .  The  following  is  an 

example  of  their  immediate  application: 


1  ^ 

J  <5?(|)  sin<|)  J  dQ  (sin(|)  cos0)  (suk})  sin6) 

•  ^  A  A 


2% 


X  exp(m  sin  (j)  cos  0  +  ib  sin  (j)  sin  0  +  ic  cos  (j))  =  . 


SPHERICALLY-ISOTROPIC  NOISE  PRESSURE  FIELD 

Let  the  incident  pressure  at  the  origin  of  coordinates,  due  to  incremental  solid  angle  d0  dcj) 
about  arrival  angle  0,(j),  be 

0,0,0)  =  exp(z(at)  P(Q,^)  t/0  d^  , 
where  noise  source  distribution  P(0,(t>)  is  zero  mean  and  satisfies 
P(0,(j))P*(0',(t>')  =  /  sin(|)  6(0  -  0')  6((|)  -  f ) . 

This  corresponds  to  an  isotropic  noise  field  of  uncorrelated  arrivals  where  /  is  a  measure  of  the 
intensity  (and  could  be  a  function  of  o),  while  the  sin(|)  factor  is  due  to  the  incremental  surface 
area  in  spherical  coordinates. 


C-2 


The  incremental  pressure  at  general  point  Pi  =  (xi,yi,zi)  is  then 
dpit,  ,  zj )  =  dp(t  -  Tj  (0,  (|)), 0,0,0) , 


where  time  delay 


Xi  sincj)  COS0  H-yj  sincj)  sm0  +  zicos 


The  total  (noise)  pressure  at  Pi  is,  therefore,  with  k  =  o/c, 

n  2% 

p(t,x„yi,z,)  =  exp(icoO  |  i/(|)  |  dQF(Q,^) 

0  0 

X  exp[/^(x,  sin  <|)  cos  0  +  y,  sin  <j)  sin  0  +  Zj  cos  <|))] . 

The  correlation  of  pressures  at  separated  times  and  distances  is 

71  271 

p(ri,x,,yi,z,)p*(?2,X2,y2,Z2)  =  exp[zro(/,  j  J0 

0  0 

It  lit 

X  J  J  dQ'  P(0,(t))P*(0',<j)')  exp[/^(x,  sin(j)  COS0  + y,  sincj)  sin0 

0  0 

+  z, cosij) -^2 sin<|)'  COS0'  — y2sin(t)'  sin0'-Z2  cos(|)')] 

7C  271 

=  I  exp[/co  (t,  -  ^2 )]  I d^  sin  <|)  J  zi0  exp[zA:({Xj  -  ^2 }  sin  (j)  cos  0 
0  0 

+{Ti  ~yi)  sin  (j)  sin  0  +  {Zj  -  Zj }  cos  <j))] 

^  T  r.  /  sin(M) 

=  47T  /  exp[z©(r,  -  ^2 )]  ^  , 

where 


d  =  [(xi  - X2)^  +  (y,  -  +  (z,  -  Z2)^]'^^ . 


This  correlation  has  the  familiar  sm{kd)l(M)  spatial  behavior,  where  d  is  the  absolute  distance 
between  the  two  points  of  interest;  this  statement  is  in  coordinate-free  form,  since  distance  d  is 
independent  of  the  particular  geometry  adopted. 
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ISOTROPIC  NOISE  VELOCITY  FIELD 


Let  w  (0,(|))  be  the  unit  vector  pointing  in  direction  0,(|): 


w^(0,<t>)  =  sin(j» COS0 T +sin<j) sin0 'f  +  cos^lc  • 

Let  the  incident  velocity  field  at  the  origin  of  coordinates,  due  to  incremental  solid  angle 
d0  d(t>  about  arrival  angle  0,<j),  be 


Jv(r, 0,0,0)  =  exp(/cor)  S(0,<|))  dQ  d^  if(0,  (j)) , 
where  noise  source  distribution  S(Q,^)  is  zero  mean  and  satisfies 
S(0,  <|))  S*  (0',  f)  =  7  sin  (j)  5(0  -  0')  5(<|)-f). 

This  corresponds  to  an  isotropic  noise  field  of  uncorrelated  arrivals  where  7  is  a  measure  of  the 
intensity  (and  could  be  a  function  of  q),  while  the  suKj)  factor  is  due  to  the  incremental  surface 
area  in  sphericail  coordinates. 

The  incremental  velocity  at  general  point  Pi  =  {xi,yi^i)  is  then 


dy(t,  xi,yi,zi)=  dy(t-Xi(e,  (j)), 0,0,0) , 
where  the  time  delay  is 

T|(0,(j))  =  (xi  sin<j)  CosO  +  yi  sm(j)  sinO  +  Z]  cos 
c  V 

The  total  (noise)  velocity  at  Pi  is,  therefore,  with  k  =  co/c. 


7C  2n 

^(t,x^,yi,z^)^exp(i(i>t)  |  J  c®S(0,(j)) 

0  0 

X  exp[/A:(xi  sin  (j)  cos  0  +  j;,  sin  ([)  sin  0  +  Zj  cos  (j))]  w(0,  ([►) . 

If  exp[  ]  denotes  the  spatial  exponential  term  above,  and  the  velocity  vector  is  expanded 
as  V  =  T  +  Vy  ^  +  Vj  ^ » then  the  velocity  components  at  Pi  are 
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Tt  In 


. (t, Xi,yi,zi)  =  exp(mt)  j*  i/cj)  |  c/0  S(0, (}>) exp[ ] sin (j)  cos 0 


0  0 


n  2n 


,  (t,  Xj ,  >>1 ,  zj )  =  exp(z(or)  J  £/(|)  J  f;?0  S(0,  (j))  exp[  ]  sin  (j)  sin  0 , 


0  0 


71  2n 


, (/, Xx,yi,z{)  =  exp(/ci)/)  J  </())  j*  J0  S(0, (j))  exp[ ] cos <j) . 


0  0 


We  are  interested  in  all  six  possible  correlations  at  general  separated  times  and  distances,  namely 


VV  vv  vv  vv  vv 

^ z  ?  ^ z  5  ^ x^ X  5  y  y  ’  ^ z^ z  ? 

where  an  abbreviated  notation  has  been  employed  in  the  latter  line. 


The  cx^^-correlation  is 


_ _  TC  Z7t 

V;t(^i>^i>Ti»^iK(^2’^2>T2>^2)  =  ®xp[zco(t,  -/j)]  J  ^/<t)  J  dQJsm^ 


X  (sin(|)  COS0)  (sin(j)  sin0)  exp[zA:({x,  -Xj}  suKj)  cos0 
+  {Ti  -  T2  }  sin  (j)  sin  0  +  {z^  -  z^ }  cos  (|))]  =  -4tz  I  exp[z(o(t]  - )] 


=  47r  /  exp[z(o(t,  -  )] 


{Q^  -3)smQ  +  3QcosQ 


k^x^-x^Xy^-y^) 


where 


=  -471  /exp[/co(t,  - 1^)]  k^{x,  -x^Xy,-y^), 


Q  =  kd  =  k  (x^-  X2f  +(yi  -  y2f  +  (z,  -  Z2f 


Here,  the  delta  function  properties  of  the  correlation  function  of  random  jSeld  S(0,(j))  and  the 
double  integral  results  derived  above  are  used.  The  ratio 72(0/0^  is  1/15  at  g  =  0. 


C-5 


It  should  be  observed  that  the  xy  correlation  involves  the  factor  (xi  -  X2)  (yi  -  y2),  in 
addition  to  the  symmetric  function  Q  of  all  the  spatial  differences.  Therefore,  by  symmetry,  the 
xz  correlation  is  identical  to  that  above,  except  that  factor  (xi  -  X2)  (y\  -  >2)  is  replaced  by 
(xi  -  X2)  (71  -  72).  A  similar  replacement  of  this  factor  holds  for  the  yz  correlation. 


The  other  type  of  (auto)  correlation  of  interest  can  be  represented  by  the  quantity 

_  71  2% 

v,(^„:ti,yi,2i)v*(/2,^2»>'2>'2’2)  =  exp[/(0(r,  -tj)]  j  cf<|)  J  ^^0  /sincj) 


X  (sin (j)  cos 9)^  exp[/^({X]  -  Xj }  sin (}>  cos 0  +  {y,  -  >^2 }  ^i^ <])  sin 0  +  {Zj  -  z^} cos (|»)] 

=  -47iF^^ /exp[?©(r, -^2)] 


=  E 


=  E 


sin(2  cosQ  (0^  - 3) sin 0  +  30 cos g  ,  2/„  „  s2 


V  0^ 


k  (x,  -X2)' 


where 


E  =  4n  I  exp 

Again,  by  use  of  the  symmetric  behavior  of  0,  the  yy  correlation  follows  immediately  upon 
replacement  of  the  (X|-X2)^  factor  by  (yi-yo)^.  A  similar  procedure  holds  for  Z7  correlation. 

When  points  Pi  and  P2  coincide,  0  =  0  and  all  three  (auto)  correlations  reduce  to  E/3; 
here,yj  (0)/0  =  1/3  is  used  at  0  =  0.  This  observation  is  consistent  vdth  the  fact  that  the  field  is 
isotropic  and  no  coordinate  direction  is  preferred,  in  so  far  as  one  point  in  space  is  concerned. 


ONE  POINT  IN  SPACE 

Now,  some  of  the  general  results  above  will  be  specialized.  For  one  point  in  space,  set  X2 
=  Xhy2=yh^2=Zh  getting  0  =  0  and 
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V;.  (?„  ^1 ,  J'l,  )v;  (rj ,  :v, ,  )  =  0 


V,  (^1 ,  ,  z, )  v;  (^2 ,  X, ,  j;, ,  z, )  =  0 

V  j,  (^1 ,  JJ'i ,  ^1 K  (?2 ,  ,  :);i ,  z, )  =  0 


\  for  all  t^,t^,Xy,y^,Zy 


Also,  there  follows 


E  An 


^w(.h^X\^yi,zi)\^{t2,x^,yx,zf^  =  —  =  — I exp[i(o (^i  - ^2 )] , 


where  subscript  w  =  x  or  3;  or  z. 

Since  the  coordinate  axes  are  arbitrary,  these  properties  can  be  summarized  in  a 
coordinate-free  form:  for  any  orthogonal  cartesian  coordinate  system,  the  three  velocity 
components  at  one  spatial  point  are  uncorrelated  with  each  other  for  all  time  delays.  Also,  each 
velocity  component  has  identical  time  (auto)  correlations. 


TWO  DIFFERENT  POINTS  IN  SPACE 


Consider  two  points  Pi  and  P2  in  space.  If  an  arbitrary  cartesian  coordinate  system  is 
erected,  the  two  locations  Pi  =  (xi,yi,zi)  and  P2  =  (X2,y2^2)  will  not  have  any  special  properties; 
that  is,  X2  5^  ^1jP2  ^  yh  ^2^  zi.  Then,  none  of  the  correlations  above  (auto  or  cross)  will  be 
zero.  Physically,  this  says  that,  generally,  all  of  the  three  orthogonal  velocity  components  at  Pi 
are  correlated  with  all  of  the  three  orthogonal  velocity  components  at  Pj. 


However,  there  is  one  special  coordinate  choice  which  does  have  uncorrelated  velocity 
components.  Since  the  noise  field  is  isotropic,  there  is  no  loss  of  generality  in  talcing  the  z-axis  to 
be  parallel  to  the  line  joining  points  Pi  and  P2;  this  still  leaves  rotational  freedom  for  the 
locations  of  the  perpendicular  x-  and  y-axes.  Thus,  y-ft\,x\,yi^i)  is  called  the  in-line  velocity 
component  at  point  Pi  at  time  ?i,  i.e.,  the  velocity  at  Pi  along  the  line  joining  points  Pi  and  P2. 
Then,  Vx(?i,xi,yi,zi)  and  Vy(rip:i  j;i,zi)  are  called  the  two  in-plane  velocity  components  at  Pi  at 
time  t\,  namely,  the  orthogonal  velocities  in  the  plane  perpendicular  to  the  z-axis  at  location 
z  =  zi .  Recall  that  the  location  of  the  x-  and  y-axes  is  arbitrary. 


Now, 

^2  =  .  y2  =  y\’ 


Q  =  k\zx -Z2\  =  kd, 


where  d  is  the  distance  between  point  Pi  and  P2.  There  follows 
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V;,  (A ,  JCi ,  3^1 ,  )v;  (^2 ,  x„  ,  Zj )  =  0 


=  Q  ^  for  all  t^,t^,x^,y^,Zy  Z. 


z,  K  (^2 ,  Xp  :Vi ,  Zj )  =  0 


Stated  in  a  coordinate-free  form,  the  in-line  velocity  component  at  P2  is  uncorrelated  with  ^  in¬ 
plane  velocity  components  at  Pi  for  all  time  delays  and  distance  separations  d.  Also,  any  in¬ 
plane  velocity  component  at  P2  is  uncorrelated  with  its  orthogonal  in-plane  velocity  component 
at  Pi  for  all  time  delays  and  distance  separations  d. 

In  addition,  for  all  , t^, x^,y^,Zy z^,  the  correlations  are 


sing  cosg')  „7i(g) 


=  E 


Q  ’ 


=  P 


sing  cosg^  _  7i(g) 

V  g^  )  Q  ’ 

-  2)  sin  g  +  2g  cos  g'* 


7,(g)  . 


V  g 


-72(0  =^Uo(0-2 


jm 

Q 


where  Q  —  kd  and  E  =  An  I  exp[2to(ri-r2)].  That  is,  the  orthogonal  in-plane  velocity  components 
have  identical  spatial  correlations,  whereas  that  for  the  in-line  velocity  component  is  different. 

All  three  (auto)  correlations  reduce  to  E/3  =  4n/3  /exp[i(D(ri-r2)]  at  Z2  =  zi,  making 
P2  =  Pi,  a  common  point.  On  the  other  hand,  for  large  physical  separations  d,  the  in-plane 
correlations  behave  as  -E  cosg/g^,  whereas  the  in-line  correlation  decays  slower,  namely,  as 
E  sing/g.  This  latter  form  is  identical  to  the  pressure  correlation  in  an  isotropic  field. 


C-8 


REFERENCES 


1 .  G.  L.  D'Spain  et  al.,  “Initial  Analysis  of  the  Data  from  the  Vertical  DIFAR  Array,” 
Proceedings  of  Mastering  the  Oceans  Through  Technology,  (OCEANS  92),  Newport,  Rhode 
Island,  October  26-29,  1992. 

2.  A.  Nehorai  and  E.  Paldi,  "Acoustic  Vector-Sensor  Array  Processing,"  IEEE  Trans,  on  Signal 
Processing,  vol.  42(9),  September  1994. 

3.  M.  Hawkes  and  A.  Nehorai,  "Surface-Mounted  Acoustic  Vector-Sensor  Array  Processing," 
IEEE  International  Conference  on  Acoustics,  Speech  and  Signal  Processing  (ICASSP96), 
May  1996. 

4.  R.  Kneipfer,  "Array  Processing  for  Detection  of  Narrowband,  Planewave  Velocity  Fields," 
Unpublished  Memorandum,  Naval  Undersea  Warfare  Center  Detachment,  New  London,  CT, 
March  1996. 

5.  R.  F.  Keltic,  “Response  of  Elastically-Coasted  Plates  to  Incident  Acoustic  Plane  Waves,” 
IPA  Activity  Report,  North  Carolina  State  University,  Raleigh,  NC,  1  November  1996. 

6.  A.  H.  Nuttall  and  B.  A.  Cray,  "Efficient  Calculation  of  Directivity  Indices  for  Certain  Three- 
Dimensional  Arrays,”  NUWC-NPT  Technical  Report  1 1,129,  Naval  Undersea  Warfare 
Center  Detachment,  New  London,  CT,  26  July  1996. 

7.  1.  S.  Gradshte3m  and  I.  M.  Ryzhik,  Table  of  Integrals,  Series,  and  Products,  Academic  Press, 
Inc.,  New  York,  NY,  1980. 


R-l/R-2 
Reverse  Blank 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

Naval  Sea  Systems  Command  (PEO-USW-ASTO,  CDR  J.  Polcari,  M.  Traweek)  2 

Office  of  Naval  Research  (334,  V.  Simmons;  331SS,  R.  Varley,  K.  Dial; 

321  US,  J.  Tague)  4 

Bolt,  Beranek  and  Newman  (N.  Martin,  R.  Brown)  2 

Litton  Guidance  and  Control  Systems  (S.  Martin)  1 

North  Carolina  State  University  (R.  Keltie)  1 

University  of  Illinois-Chicago  (A.  Nehorai)  1 

Scripps  Institution  of  Oceanography  (G.  D’Spain)  1 

Defense  Technical  Information  Center 


2 


